Sugar-sweetened beverages, low/no-calorie beverages, fruit juice and non-alcoholic fatty liver disease defined by fatty liver index: the SWEET project

Background Sweetened beverage intake may play a role in non-alcoholic fatty liver disease (NAFLD) development, but scientific evidence on their role is limited. This study examined associations between sugar-sweetened beverages (SSB), low/no-calorie beverages (LNCB) and fruit juice (FJ) intakes and NAFLD in four European studies. Methods Data for 42,024 participants of Lifelines Cohort, NQPlus, PREDIMED-Plus and Alpha Omega Cohort were cross-sectionally analysed. NAFLD was assessed using Fatty Liver Index (FLI) (≥60). Restricted cubic spline analyses were used to visualize dose–response associations in Lifelines Cohort. Cox proportional hazard regression analyses with robust variance were performed for associations in individual cohorts; data were pooled using random effects meta-analysis. Models were adjusted for demographic, lifestyle, and other dietary factors. Results Each additional serving of SSB per day was associated with a 7% higher FLI-defined NAFLD prevalence (95%CI 1.03–1.11). For LNCB, restricted cubic spline analysis showed a nonlinear association with FLI-defined NAFLD, with the association getting stronger when consuming ≤1 serving/day and levelling off at higher intake levels. Pooled Cox analysis showed that intake of >2 LNCB servings/week was positively associated with FLI-defined NAFLD (PR 1.38, 95% CI 1.15–1.61; reference: non-consumers). An inverse association was observed for FJ intake of ≤2 servings/week (PR 0.92, 95% CI: 0.88–0.97; reference: non-consumers), but not at higher intake levels. Theoretical replacement of SSB with FJ showed no significant association with FLI-defined NAFLD prevalence (PR 0.97, 95% CI 0.95–1.00), whereas an adverse association was observed when SSB was replaced with LNCB (PR 1.12, 95% CI 1.03–1.21). Conclusions Pooling results of this study showed that SSB and LNCB were positively associated with FLI-defined NAFLD prevalence. Theoretical replacement of SSB with LNCB was associated with higher FLI-defined NAFLD prevalence. An inverse association was observed between moderate intake of FJ and FLI-defined NAFLD. Our results should be interpreted with caution as reverse causality cannot be ruled out.


INTRODUCTION
Non-alcoholic fatty liver disease (NAFLD) is a general term that refers to a broad range of liver disorders, including benign macrovesicular hepatic steatosis (>5% of hepatocytes), nonalcoholic steatohepatitis, hepatic fibrosis, and liver cirrhosis, which constitutes the third most common cause of hepatocellular carcinoma [1]. Globally, NAFLD prevalence is estimated to be around 25%, which has steadily increased and currently coincides with obesity, type 2 diabetes, and metabolic syndrome epidemic [2].
Major determinants of NAFLD include sedentary lifestyles and poor-quality diets, namely high sugar intake [3][4][5]. A recent metaanalysis of 12 studies (nine cross-sectional, two case-control and one cohort study) showed that sugar-sweetened beverages (SSB) intake was associated with a 39% (95% CI 1.29-1.50) higher NAFLD risk [6]. However, studies on the association between SSB intake and NAFLD were mainly conducted in relatively small studies using liver imaging techniques or biopsies to define NAFLD, which are often not feasible in large population studies [6,7]. Surrogate formulas-employing serum biomarkers-have been developed, validated and widely used in recent decades, including the Fatty Liver Index (FLI) by Bedogni et al. [8][9][10], to facilitate non-invasive and simple NAFLD assessment in larger studies. However, evidence on the association between SSB and NAFLD detected with FLI in larger studies is still limited [3,11,12].
Low/no-calorie sweeteners have been increasingly used as an alternative to sugar. However, evidence on the association between low/no calorie beverages (LNCB) and NAFLD is limited and generally shows no association [7,13,14]. Besides LNCB, fruit juice (FJ) is also often considered an alternative to SSB, although some argue that FJ contains high amounts of sugar apart from its vitamin and polyphenol content. Population-based studies showed beneficial associations with cardiometabolic health at moderate doses of FJ consumption [15][16][17], but evidence specifically on NAFLD is lacking.
Given all the aforementioned considerations, we aimed to investigate the association between SSB, LNCB, and FJ intakes and FLI-defined NAFLD in four European studies. Additionally, we investigated whether LNCB and FJ could be used as an alternative to SSB by means of substitution analyses which is currently lacking in the literature.

Study design and population
The SWEET project is a European Union-funded project that aims to investigate public health and safety, obesity, and sustainability risks and benefits of replacing sugar with sweeteners and sweetness enhancers (www.sweetproject.eu). The results presented in this manuscript are based on harmonized data of the Lifelines Cohort study (Lifelines; The Netherlands), the Nutrition Questionnaire Plus study (NQPlus; The Netherlands), the PREDIMED-Plus study (PREDIMED-Plus; Spain), and the Alpha Omega Cohort (AOC; The Netherlands). All participants in all cohorts gave written informed consent before participating.

Lifelines
Lifelines is a multidisciplinary prospective population-based cohort study examining in a unique three-generation design the health and healthrelated behaviours of 167,729 persons living in the North of The Netherlands [18,19]. It employs a broad range of investigative procedures in assessing biomedical, sociodemographic, behavioural, physical and psychological factors that contribute to the health and disease of the general population, with a special focus on multi-morbidity and complex genetics. Between 2006 and 2013, participants aged 0-93 yr were recruited to undergo baseline measurements. For present analysis, data for 152,728 participants aged ≥18 yr were included. After consecutive exclusion of those with missing dietary data (n = 8633), implausible energy intake (<800 or >4000 kcal/day for men or <500 or >3500 kcal/day for women; n = 15,483) [20], missing FLI data (n = 86,137), history of hepatitis (n = 31), excessive alcohol consumption (>20 g/day for women and >30 g/day for men; n = 1825), or missing covariates (n = 3679), n = 36,940 remained for current analyses (Supplemental Table 1). Lifelines has been approved by the Medical Ethical Review Committee of the University Medical Center in Groningen.

NQPlus
NQPlus is a prospective cohort study involving mainly Caucasian Dutch adults aged 20-70 yr from the central part of the Netherlands (Wageningen, Ede, Renkum, Arnhem, Barneveld, Veenendaal) [21]. NQPlus' objectives were to establish a database for national dietary assessment reference to produce and validate food frequency questionnaires (FFQs), as well as to provide a database for longitudinal study in dietary factors and healthrelated outcomes. Extensive data collection was performed on habitual dietary intake and various health outcomes, including cardiometabolic health parameters. Between 2011 and 2013, 2048 participants were included in the study. After participants with missing dietary data (n = 401), implausible energy intake (n = 20), missing FLI data (n = 56), history of hepatitis (n = 32) or excessive alcohol consumption (n = 184) were consecutively excluded, n = 1341 remained (Supplemental Table 1). Due to the high percentages of missing covariates (n physical activity = 110 [8%], n sedentary behaviour = 111 [8%], n smoking status = 126 [9%]), multiple imputations were applied using "mice" package in R and five datasets were produced [22]. NQPlus was approved by the ethical committee of Wageningen University and Research.

PREDIMED-Plus
PREDIMED-Plus is a multicentre randomized controlled trial which was conducted based on the idea of the previous study, the PREDIMED (in Spanish: PREvención con DIeta MEDiterránea) trial [23]. PREDIMED-Plus aims to examine the effect of a traditional Mediterranean diet with energy restriction and increased physical activity on weight reduction and cardiovascular incidence and mortality. PREDIMED-Plus was conducted in 23 study centres across Spain. From September 2013 to November 2016, PREDIMED-plus recruited men aged 55-75 yr and women aged 60-75 yr with overweight or obesity who met at least three metabolic syndrome criteria according to the International Diabetes Federation and the AHA/ National Heart, Lung and Blood Institute [24]. For present analysis, data from the Canary Island study centre were used, which included information for 266 participants with complete dietary data. Participants with implausible energy intake (n = 3), missing FLI (n = 10), and excessive alcohol consumption (n = 1) were consecutively excluded, and n = 250 remained for the current analyses (Supplemental Table 1). PREDIMED-Plus was registered at the International Standard Randomized Controlled Trial database (ISRCTN; 89898870) and the protocol was approved by local Research Ethics Committees.

AOC
AOC is a prospective cohort study of 4837 Dutch patients aged 60-80 yr who had experienced a myocardial infarction (MI) within 10 yr before study enrolment [25,26]. The cohort is continuously followed for cause-specific mortality. Baseline measurements took place in 2002 and 2006 during which patients filled in questionnaires on demographics, lifestyle, diet, medication use and medical history, and underwent physical examination by trained research nurses. For current analyses, patients with missing dietary data (n = 453) or implausible dietary intake (n = 24), missing FLI (n = 179), excessive alcohol consumption (n = 647), or missing covariates (n = 41) were excluded, resulting in a total of n = 3493 remaining for analysis (Supplemental Table 1). AOC was registered at ClinicalTrials.gov (NCT03192410) and was approved by a central Medical Ethics Committee in the Netherlands.

Dietary assessment
Dietary intake data were collected using a newly developed and validated 110-item FFQ [27] and 183-item semi-quantitative FFQ [21,28] in Lifelines and NQPlus, respectively, which covered all major food groups. These FFQs were created based on the Dutch consumption survey and were validated for energy and a range of nutrients as well as intake biomarkers. Average daily nutrient intake was calculated by multiplying consumption intake frequency by portion size and nutrient content in grams as indicated in the Dutch Food Composition Table (NEVO) 2011. In PREDIMED-Plus, dietary intake was collected using a 143-item FFQ [23,29,30], based on usual food consumption in Spain and designed to assess overall diet. This FFQ has been validated in different studies using four 3-or 4-day food records. Average daily nutrient intake was calculated by multiplying consumption intake frequency by portion size and nutrient content in the Spanish food composition tables. Within AOC, dietary data were collected using a semiquantitative 203-item FFQ, a modified and extended version of the validated FFQ that was specifically designed to estimate fatty acids and cholesterol intake [31]. The Dutch Food Composition (NEVO) 2006 was used to convert food intake to energy and nutrient intakes in AOC. All population studies include intake measures of SSB, LNCB, and FJ, except for AOC where LNCB was combined with water, which therefore could not be analysed individually. Definition of SSB and LNCB covers soft drinks or lemonade sweetened with sugar or sugar substitutes. Coffee or tea sweetened with sugar or sweeteners as well as sweetened dairy drinks were not included in the definition of SSB and LNCB. FJ was defined as fruit juice i.e., apple juice, orange juice and mixed fruit juice, which were mainly pasteurized. SSB, LNCB, and FJ intakes were presented as servings of 150 ml. well as liver enzymes including gamma-glutamyl transferase (GGT) were measured using routine procedures on a Roche Modular P800 chemistry analyser in Lifelines, Dimension Vista 1500 automated analyser or a Roche Modular P800 chemistry analyser in NQplus, and Cobas 8000 c 702 Module Roche Diagnostics in PREDIMED-Plus. In AOC, TG was analysed using Roche Hitachi 912, whereas liver enzymes were determined using automated analyses Abbott Architect ci8200. Anthropometry measurements (weight, height and waist circumference) were performed by well-trained staff. Body weight and height were measured without shoes. Body mass index (BMI) was calculated as weight divided by squared height (kg/m 2 ). NAFLD was defined using FLI ≥ 60, which has been established by Bedogni et al. (2006) [8]. FLI is calculated based on serum concentrations of TG and GGT, BMI, and waist circumference, applying the following formula:

Assessment of covariates
Information on sociodemographic, lifestyle, and disease history was obtained from self-or interviewer-administered questionnaires. Educational levels was categorized into low (less than secondary qualification), medium (from secondary qualification up to university with no bachelor's degree) or high (Bachelor's, Master's, or Ph.D. degree). Smoking status was categorized as non-smoker, former smoker, or current smoker. In Lifelines and NQplus, physical activity and sedentary behaviours (TV-watching) was assessed using the Short Questionnaire to Assess Health (SQUASH) [32] and the Activity Questionnaire for Adults and Adolescents (AQuAA) [33], respectively. In PREDIMED-Plus, physical activity was assessed using the Physical Activity Readiness Questionnaire (PAR-Q) as well as the Rapid Assessment of Physical Activity questionnaires 1 and 2 (RAPA-1 and RAPA-2) [34]. Sedentary activity in PREDIMED-Plus was assessed by using the Nurses' Health Study physical activity questionnaire validated for the Spanish population [35]. Physical activity in Lifelines, NQPlus and PREDIMED-Plus was reported in Metabolic equivalent (MET)-min/week for moderate-level activity and in min/week for sedentary behaviour. In AOC, physical activity was assessed by using validated the Physical Activity Scale for the Elderly [36], and categorized as low (no activity/light activity, ≤3 METs), moderate (0-5 day/week moderately or vigorously active, >3 METs) or high (≥5 day/week moderately or vigorously active, >3 METs) physical activity. Alcohol intake was assessed using FFQ from which ethanol intake was calculated and categorized as 0, >0-≤10, >10-≤20, or >20 g/day.

Statistical analysis
Baseline characteristics were presented as mean with standard deviation for normally distributed continuous variables or as a median and interquartile range for variables with skewed distribution. Categorical variables were shown as numbers and percentages. Dose-response associations between SSB, LNCB, and FJ intakes and FLI-defined NAFLD prevalence were first analysed using restricted cubic spline analyses (3 knots) [37]. To ensure adequate power and precision, restricted cubic spline analyses were only performed for Lifelines. The fit of the spline model was examined against a linear model with the likelihood-ratio test. Cox proportional hazard regression with robust variance estimate was used to investigate the associations of SSB, LNCB and FJ intakes with FLI-defined NAFLD prevalence with a 95% confidence interval (PR [95% CI]) in each cohort, whereas theoretical substitution analyses were conducted using the leave-one-out model, where the model included SSB, LNCB, and FJ (in servings/day) as one variable followed by beverage defined as a replacement. Results of individual cohorts were subsequently pooled using random effect meta-analyses. To explore the potential of reverse causality, sensitivity analyses were conducted in Lifelines after excluding those with self-reported medical conditions (i.e., diabetes, CVD, hypercholesterolemia, or hypertension). All analyses were adjusted for age, sex (Model 1), educational levels (low, medium, or high), lifestyle factors including moderate physical activity (MET-min/week or in categories as in AOC), sedentary behaviour (min/week), smoking status (never, former, or current smoker), and alcohol use (0, >0-≤10, >10-≤20, or >20 g/day) (Model 2), intakes of other food groups including grains (g/day), potatoes (g/day), vegetables (g/day), fruit (g/day), meat and processed meat (g/day), dairy (g/d), coffee (ml/day), tea (ml/day), legumes (g/day), nuts (g/day), oils and fats (g/day), sugary foods (candy, sweet snack i.e., cookies, milk chocolate, any kind of toppings and popsicle) (g/day), and mutual adjustment for other beverages (SSB, LNCB, and FJ in g/day) as well as total energy intake (kcal/d) (Model 3). Statistical analysis was performed using R 4.0.2 and RStudio 1.3.959 for NQplus, PREDIMED-Plus, and AOC, and RStudio 2022.02.0 for Lifelines.
In most cohorts, participants with a higher SSB intake were more likely to be younger, men and less physically active. The opposite trends were observed in AOC. In all cohorts, SSB consumers were more likely to be diagnosed with diabetes, had a higher intake of meat, sugary food and total energy, and lower intake of vegetables and fruit (Supplemental Table 2). Participants with higher LNCB intake were more likely to be younger, less physically active, to be diagnosed with diabetes, hypercholesterolemia or hypertension, have a BMI ≥25, higher meat intake, and lower intake of vegetables and fruit (Supplemental Table 3). For FJ, participants with a higher intake were more likely to have BMI <25 and higher sugary food consumption (Supplemental Table 4).

Association between SSB, LNCB and FJ intakes and FLIdefined NAFLD prevalence
Dose-response analysis in Lifelines did not provide strong evidence of a non-linear association between SSB and FLIdefined NAFLD (P non-linearity = 0.07) (Fig. 1). Pooled results of four studies showed an association between each additional serving/day of SSB and higher prevalence of FLI-defined NAFLD (PR 1.07, 95% CI 1.03-1.11) ( Table 2). Although there was evidence of interaction between SSB and age categories, results from a stratified analysis in Lifelines by age as well as by sex did not indicate significant differences (Supplemental Table 6). Sensitivity analysis by excluding those with self-reported disease in Lifelines did not alter the associations (Supplemental Table 7).
Dose-response analysis in Lifelines showed strong evidence of a non-linear association between LNCB and FLI-defined NAFLD prevalence (P non-linearity = <0.001) (Fig. 1), with the association getting stronger when consuming ≤1 serving/day and levelling off at higher intake levels. Subsequent Cox regression showed that, compared to no intake, LNCB intake of ≤2 servings/week was associated with a PR of 1.13 (95% CI 1.00-1.26), whereas intake of >2 servings/week was associated with a PR of 1.38 (95% CI 1.15-1.61) ( Table 3). In Lifelines, stratified analysis by sex showed a stronger association in women (PR 1.65, 95% CI 1.54-1.77) than in men (PR 1.36, 95% CI 1.29-1.43) when comparing those with LNCB intake >2 servings/week with non-consumers (Supplemental Table  6). Excluding participants with self-reported disease did not alter the main findings (Supplemental Table 7). Theoretical replacement of SSB with the same amount of LNCB was associated with a higher prevalence of FLI-defined NAFLD (PR 1.12, 95% CI 1.03-1.21) ( Table 4).
Dose-response analysis in Lifelines did not suggest evidence of a nonlinear association between FJ and FLI-defined NAFLD prevalence (P non-linearity = 0.86) (Fig. 1). However, pooled analysis by categories showed an inverse association between FJ and FLIdefined NAFLD at intake levels of ≤2 servings/week (PR 0.92, 95% CI 0.88-0.97), but no association was seen at higher intake levels (PR 1.00, 95% CI 0.97-1.04) when compared to no intake (Table 3). After excluding those with self-reported disease, the associations did not change significantly (Supplemental Table 6). Finally, there was no

DISCUSSION
Harmonized data analyses of four European studies encompassing both general population and cardiometabolic subjects showed positive associations between SSB and LNCB intakes and FLI-defined NAFLD prevalence. Dose-response analyses in Lifelines showed a non-linear positive association between LNCB and FLI-defined NAFLD with the stronger association at intake levels of <1 serving/day and levelling off at higher intake levels.
A beneficial association was observed for moderate FJ intake ranging from ≤2 servings/day when compared to no intake. Theoretical replacement of SSB with the same amount of LNCB showed a positive association with FLI-defined NAFLD, whereas no significant association was observed when replacing SSB with FJ. Overall, our findings on the association between SSB and FLIdefined NAFLD prevalence are consistent with previous research [6,7,11,38,39]. A meta-analysis of 12 studies (nine crosssectional, two case-control and one cohort study) showed a 39% (95% CI 1.29-1.50) higher NAFLD risk-as assessed using liver imaging or biopsy-when comparing highest SSB intake group with non-consumer [6]. Furthermore, a recent study by Zhang et al. in 14,845 participants showed a 47% (HR 95% CI 1.25-1.23) and 59% (HR 95% CI 1.07-2.37) higher prevalence of NAFLD when NAFLD was assessed using abdominal ultrasonography and serum biomarkers (hepatic steatosis index) [11]. In addition, results from short-term trials, experimental studies, and meta-analyses of RCT evaluating individual metabolic markers further confirm our findings [40][41][42][43][44]. The biological explanation of the adverse association between SSB intake and NAFLD can be explained by  Fig. 1 Dose-response associations of SSB, LNCB and FJ intakes with FLI-defined NAFLD in Lifelines. Three knots with 0 g/day as a reference value were placed and model was adjusted for age, sex, education levels, moderate physical activity, sedentary behaviour, smoking status, alcohol intake, grains, potatoes, vegetable, fruit, meat and processed meat, dairy, coffee, tea, legumes, nuts, oils and fats, sugary foods, mutual adjustment for other beverages (SSB, LNCB, or FJ), and energy intake. CI confidence interval, FJ fruit juices, LNCB low/no-calorie sweetened beverages, PR prevalence ratio, SSB sugar-sweetened beverages.  (27) 2897 (20) 1454 (18) 3849 (27) 2180 (27) 2322 (20) 3698 (  several biological mechanisms. Its high fructose content may induce hepatic de novo lipogenesis and possibly insulin resistance [45,46]. Moreover, lack of energy compensation after liquid calorie intake affects energy balance disruption, which contributes to weight gain, an important factor in NAFLD development [47]. We observed positive associations between LNCB, as well as between theoretical replacement of SSB with LNCB and FLI-defined NAFLD prevalence. To date, evidence on LNCB and NAFLD is still limited [7,13,14], especially from large studies. In line with our findings, prospective analysis among 1636 women participating in the Framingham Heart Study showed a 48% higher risk of NAFLD and LNCB-when NAFLD was detected using computed liver attenuation measurements [7]. However, this association disappeared after further adjustment for BMI (OR 1.11, 95% CI 0.80-1.56). A prospective study in women with a history of gestational diabetes also showed no association between LNCB and the majority of NAFLD markers when the models were adjusted by various covariates, particularly pre-pregnancy BMI, as well as when performing sensitivity analyses excluding those with chronic diseases -indicating that reverse causality may be the cause [14]. Within our study, participants at higher intake levels of LNCB tended to have a BMI ≥ 25, which might also indicate the presence of reverse causality. Participants with NAFLD or elevated NAFLD parameters may have switched from SSB to LNCB to control their health, e.g., weight regulation. Evidence from experimental studies was summarized in a meta-analysis of 12 RCTs involving 601 adults with overweight or obesity showing that replacing SSB with LNCB reduced body weight, body fat, BMI, intrahepatocellular lipid, but not other clinical parameters including liver enzymes i.e. alanine transaminase and aspartate aminotransferase [41]. A biological explanation for the association between LNCB with health outcomes is inconclusive. LNCB has been suggested to disturb reward systems, activate the cephalic phase insulin response, and induce gut microbiota dysbiosis, which can further lead to insulin resistance and the development of type 2 diabetes [48,49]. However, more studies including experimental studies in humans as well as long-term evaluations in large cohorts are still needed to further investigate these hypotheses [7,50].
To date, only a few studies specifically explored associations of FJ with health outcomes. Our results suggest an inverse association when consuming ≤2 servings/day compared to no intake, which has also been observed in a dose-response meta-analysis of two cohorts investigating the association between FJ intake (75-150 ml/day) and MetS incidence [15]. Khan et al. also demonstrated an inverse association between FJ intake of~150 ml and CVD incidence, but not at higher intake levels [51]. A recent and extensive meta-analysis by D'Elia also suggested a non-linear association between low to moderate 100% FJ intake (<80 ml/day) with incidence rates of stroke and CVD [17]. This inverse association among participants with moderate intake levels might be due to overall healthier behaviours of participants within this intake group. However, this notion is not clearly seen in our study populations as the distribution of characteristics was comparable across groups in all studies. Possible underlying mechanisms of the inverse associations between moderate intake of FJ and NAFLD include the beneficial composition of FJ, which may contain a certain amount of antioxidants (i.e polyphenols) and other bioactive components (i.e vitamins, fibre and minerals) that enhance metabolic profile [15,17,51,52]. However, its high caloric content and the presence of fructose may be counterproductive to the benefits of consuming higher quantities of FJ.
Strengths of our study include the large sample size derived from several European population-based studies, the opportunity to control for a wide range of relevant confounders and to perform stratified analysis particularly in Lifelines which comprised the largest cohort. We conducted the harmonized statistical analyses across cohorts including substitution analyses which to date, are often lacking in recent studies. This is the first metaanalysis focusing on NAFLD assessed by FLI, which facilitated noninvasive and simple assessment in larger studies. There are also several limitations to this study. Dietary intake was self-reported, hence bias due to social desirability cannot be excluded. Although the FFQs used in current study were extensively evaluated for a range of nutrients and foods, they were not specifically designed to examine SSB, LNCB, and specific sweeteners intakes. Moreover, SSB, LNCB, and FJ intakes in our study were relatively low, which may explain the lower risk estimates in our study when compared to some previous ones i.e from the US [6,7,13]. SSB intake in our Dutch study populations was lower than that reported in the Dutch National Food Consumption Survey, which may be because SSB was grouped together with carbonated/soft/isotonic drinks and diluted syrups. We used FLI as a proxy for NAFLD, which has limitations including its inability to distinguish moderate to severe steatosis from mild steatosis [9,53]. The gold standard for NAFLD diagnosis is liver biopsy which is an invasive and costly procedure.
To avoid these consequences, non-invasive imaging procedures have been introduced, i.e ultrasound or magnetic resonance. However, it also has limitations associated with availability and cost for large-scale studies, therefore use of serum biomarkers is preferred. FLI has been externally validated in several European and Asian populations [10,[54][55][56][57], and has been endorsed to be the best-validated tool to detect the presence of steatosis [58,59]. Despite adjustment for a wide range of confounders, residual confounding may still be present. However, additional sensitivity analysis using the E-value [60] showed that any unmeasured confounder would need to be quite strongly associated with both exposure and outcome to explain away the association (data not shown), which is not very likely based on the results for the confounders included. Finally, as always in a cross-sectional study, reverse causality cannot be excluded.
In conclusion, we observed adverse associations between SSB and LNCB intakes and FLI-defined NAFLD prevalence, as well as between replacement of SSB with the same amount of LNCB and FLI-defined NAFLD, which may partly be explained by reverse causality. Our findings suggest a beneficial association between moderate intake of FJ and FLI-defined NAFLD at intake level of ≤2 servings/day when compared to no intake. These findings provide additional information on the potential adverse impact of sweetened beverages on health. Longer-term prospective studies with objective methods determining the intake of sugar and sweeteners are warranted to further substantiate our findings.

DATA AVAILABILITY
The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request and after approval from the relevant partners.